Method and apparatus for detecting and locating noise sources whether correlated or not

ABSTRACT

A method of detecting and locating noise sources each emitting respective signals S j  with j=1 to M, detection being provided by a sound wave or vibration sensors each delivering a respective time-varying electrical signal s i  with i in the range 1 to N, wherein the method steps include: (a) taking the time-varying electrical signals delivered by the sensors, each signal s i (t) delivered by a sensor being the sum of the signals S j  emitted by the noise sources; (b) amplifying and filtering the time-varying electrical signals as taken; (c) digitizing the electrical signals; (d) calculating a functional; and (e) minimizing the functional relative to the vectors n j  for j=1 to M so as to determine the directions of vector n j  of the noise sources.

FIELD OF THE INVENTION

The present invention relates to detecting and locating sources of noise in the general sense, using sensors that are appropriate for the nature of the noise source.

The invention relates to a method of detecting and locating noise sources disposed in a space of one, two, or three dimensions and optionally correlated with one another, and presenting emission spectra of narrow or broad band.

The invention finds particularly advantageous applications in the field of locating sources of noise optionally accompanied by echo and coming, for example, from vehicles, ships, aircraft, or firearms.

BACKGROUND OF THE INVENTION

In numerous applications, a need arises to be able to locate in relatively accurate manner a source of noise in order to take measures to neutralize it. Numerous solutions are known in the prior art for acoustically locating noise sources. The main known solutions make use of techniques for correlating signals delivered by detection sensors.

Those techniques present the drawback of being particularly sensitive to interfering noise occurring in the environment of the measurement sensors. Furthermore, it must be considered that those techniques constitute specific methods that are adapted to each application under consideration.

The technique in most widespread use involves antennas having a large number of sensors (several hundred) and a large computer system implementing beam forming so as to aim in a given direction in order to increase the signal-to-noise ratio. That method does not make any a priori assumption concerning the number of sources and any possible correlation between them, which leads to a loss of resolution.

OBJECTS AND SUMMARY OF THE INVENTION

There therefore exists a need to have a general method of detecting and locating noise sources in space, when the number of noise sources is small and is known or overestimated.

The invention seeks to satisfy this need by proposing a method of detecting and locating noise sources by means of sensors adapted to the nature of the noise source, the method presenting low implementation costs.

To achieve this object, the method of the invention consists:

-   -   in taking the time-varying electrical signals delivered by the         sensors, each signal s_(i)(t) delivered by a sensor being the         sum of the signals S_(j) emitted by the noise sources;     -   in amplifying and filtering the time-varying electrical signals         as taken;     -   in digitizing the electrical signals;     -   in calculating the functional f, such that:         ${f\left( {n_{l},\ldots\mspace{11mu},n_{j},\ldots\mspace{11mu},n_{M}} \right)} = \frac{\det\left( {{< {T_{k}(\omega)}},{{T_{l}^{*}(\omega)} > \mspace{14mu} k},{l = {0{\mspace{11mu}\;}{to}\mspace{14mu} M}}} \right)}{\det\left( {{< {T_{k}(\omega)}},{{T_{l}^{*}(\omega)} > \mspace{14mu} k},{l = {1\mspace{14mu}{to}\mspace{14mu} M}}} \right)}$         with         $\left( {T_{k}(\omega)} \right)_{i} = {\mathbb{e}}^{J\;\omega\frac{{< n_{k}},{c_{i} >}}{c}}$         -   <.,.> being the scalar product;         -   c_(i) being the vector constructed between the center of             gravity of the sensors and the position of sensor i;         -   n_(j) being the unit vector in the direction defined by the             center of gravity of the senors and source j;         -   with T₀=s; and         -   with c=the speed of sound; and     -   in minimizing the functional f relative to the vectors n_(j) for         j=1 to M in such a manner as to determine the directions n_(j)         of the noise sources.

BRIEF DESCRIPTION OF THE DRAWINGS

Various other characteristics appear from the description given below with reference to the accompanying drawing which shows embodiments and implementations of the invention as non-limiting examples.

FIG. 1 is a diagram showing the principle of the detection method of the invention.

FIG. 2 is a diagram showing a detail characteristic to the method of the invention.

MORE DETAILED DESCRIPTION

As can be seen in FIG. 1, the method of the invention consists in locating noise sources X₁, X₂, . . . , X_(j), . . . , X_(M) where j varies over the range 1 to M, the sources being distributed in space and each emitting a respective signal S_(j) with j varying in the range 1 to M. The method of the invention consists in locating the noise sources X_(j) using sound wave or vibration sensors Y₁, Y₂, . . . , Y_(i), . . . , Y_(N) where i varies over the range 1 to N, each delivering a respective time-varying electrical signal s₁, s₂, . . . , s_(i), . . . , s_(N).

The method consists in taking the time-varying electrical signals s_(i)(t) delivered by each of the sensors and representative of the sums of the signals S_(j) emitted by the noise sources X_(j). The signals s_(i)(t) received on the N senors on the basis of the sum of the contributions of the various sources is written as follows: ${s_{i}(t)} = {\sum\limits_{j = 1}^{M}\;{A_{ij}{S_{j}\left( {t - \frac{r_{ij}}{c}} \right)}}}$ where i=1 to N, r_(ij) is the distance between the noise source X_(j) and the sensor Y_(i), and c is the speed of sound in the ambient medium.

The term A_(ij) represents the attenuation due to propagation together with the sensitivity factor of the sensors and is expressed as follows: A _(ij) =B _(i) C(r _(ij)) where i=1 to N and j=1 to M, where B_(i) is the sensitivity coefficient of sensor Y_(i) and where C(r_(ij)) is the attenuation coefficient due to propagation over a distance r_(ij).

The sensors Y_(i) are associated with respective electronic units (not shown) for amplifying and lowpass filtering the signals they pick up. The sensors are preferably matched in modulus and phase so that their sensitivities are identical. Thus, B_(i)=G for i=1 to N.

Advantageously, in order to facilitate implementing the antenna of sensors as defined above, the sensors Y_(i) are placed relatively close to one another. Consequently, for remote sources, the distance r_(ij) is of the order of the distance r_(j), i.e. the distance between the center of gravity of the sensors and the source X_(j), Thus, attenuation becomes a function of the distance r_(j) only with C(r_(ij))=C(r_(j)), with i=1 to N and j=1 to M.

It can be deduced therefrom that: A _(ij) =G.C(r _(j))=a(r _(j)) where i=1 to N and j=1 to M and: ${s_{i}(t)} = {\sum\limits_{j = 1}^{M}\;{{a\left( r_{j} \right)}{S_{j}\left( {t - \frac{r_{ij}}{c}} \right)}}}$ where i=1 to N.

Since the amplitudes of the sources X_(j) are unknown, the following equation can be written as follows, integrating the term a(r_(j)) in S_(j): ${s_{i}(t)} = {\sum\limits_{j = 1}^{M}\;{S_{j}\left( {t - \frac{r_{ij}}{c}} \right)}}$ where i=1 to N.

Using Fourier transforms, the expression for the signals s_(i)(t) becomes: $\begin{matrix} {{{\hat{s}}_{i}(\omega)} = {\sum\limits_{j = 1}^{M}\;{{{\hat{S}}_{j}(\omega)} \cdot {\mathbb{e}}^{{- J}\;\omega\frac{r_{ij}}{c}}}}} & (1) \end{matrix}$ where i=1 to N where ŝ and Ŝ are the Fourier transforms of s and S respectively and where ω is angular frequency.

This first equation (1) relates the received signals to the distance r_(ij), i.e. to the positions of the sources X_(j).

As can be seen in FIG. 2, other relationships can be expressed, associated with geometrical considerations enabling the distances r_(ij) to be related to the unit vector n_(j), which determines the direction defined by the center of gravity of the sensors and the source generating the signal S_(j). The position of the senors is defined by the vector C_(i) constructed from the positions of the sensors Y_(i) and the position of their center of gravity. A development restricted to the first order of r_(ij) then provides: r _(ij) ≈r _(j) −<n _(j) ,c _(i)>  (2) where i=1 to N and j=1 to M, and where <., .> is the scalar product.

Thus, by replacing r_(ij) by the approximate expression given in (2) and integrating the phase term: ${\mathbb{e}}^{{- J}\;\omega\frac{r_{j}}{c}}$ which depends only on the source X_(j) in the magnitude Ŝ_(j)(ω), equation (1) can be written: $\begin{matrix} {{{\hat{s}}_{i}(\omega)} = {\sum\limits_{j = 1}^{M}\;{{{\hat{S}}_{j}(\omega)} \cdot {\mathbb{e}}^{{- J}\;\omega\frac{{< n_{j}},{c_{i} >}}{c}}}}} & (3) \end{matrix}$ where i=1 to N.

This relationship can also be expressed in matrix and vector form: $\begin{matrix} {{{\hat{s}}_{i}(\omega)} = {\sum\limits_{j = 1}^{M}\;{{{\hat{S}}_{j}(\omega)} \cdot {T_{j}(\omega)}}}} & (4) \end{matrix}$ with, for ith coordinate of the vector T_(j): $\left( T_{j} \right)_{i} = {\mathbb{e}}^{{- J}\;\omega\frac{{< n_{j}},{c_{i} >}}{c}}$ where i=1 to N. Or indeed: s(ω)=T.S(ω)  (5) where T=matrix having the general term: $T_{ij} = {\mathbb{e}}^{{- J}\;\omega\frac{{< n_{j}},{c_{i} >}}{c}}$

In the presence of additive noise, equation (4) becomes: $\begin{matrix} {{s_{i}(\omega)} = {{\sum\limits_{j = 1}^{M}\;{{{\hat{S}}_{j}(\omega)} \cdot {T_{j}(\omega)}}} + {B(\omega)}}} & (6) \end{matrix}$ where B is the noise vector which depends on ω.

The method of the invention consists in determining the directions of the sources X_(j) defined by the vectors n_(j) for j=1 to M.

When the sources X_(j) are arbitrary, i.e. correlated or non-correlated, the probability of the presence of Gaussian noise B at the sensors Y_(i) is given by: b ⋅ 𝕖^(−a∫B² ⋅ 𝕕ω) where a and b depend on the variance.

Thus, from equation (6), the most probable position for the source is the position which minimizes the following magnitude: ${{{s(\omega)} = {\sum\limits_{j = 1}^{M}\;{{{\hat{S}}_{j}(\omega)} \cdot {T_{j}(\omega)}}}}}^{2}$

In other words, the projection of s onto the direction orthogonal to the hyperplane generated by the vectors T_(j) for j=1 to M must be of minimum norm.

That constitutes the square of the height of the parallelepiped constructed on the vectors s and T_(j), said height h being calculated as the ratio of the volume V to the base area S, i.e.: $h = \frac{V}{S}$

The magnitudes V and S are expressed as a function of the determinants of the Gramm matrices in which the element (k, l) is constituted by the scalar product: <T_(k),T_(l)*> with T_(l)* being the vector that is the conjugate of T_(l).

Thus, S ² =det(<T _(k) ,T _(l) *>k,l=1 to M) V ² =det(<T _(k) ,T _(l) *>k,l=0 to M) in which it is assumed T₀=s.

Thus: ${B}^{2} = \frac{V^{2}}{S^{2}}$ or indeed: f(n ₁ , . . . , n _(j) , . . . , n _(M))=∥B∥ ²

This is a function of the vectors n_(j), which in three dimensions depends only on two angles θ_(j) and φ_(j), i.e. elevation and azimuth, and on angular frequency ω. Any a priori knowledge about the spectra of the sources Ŝ_(j) can also be used. For example:

-   -   for narrow band noise sources, provision is made to minimize the         following functional f₁: $\begin{matrix}         {f_{1} = {\sum\limits_{k}\;{{B(\omega)}_{k}}^{2}}} & (7)         \end{matrix}$     -    ω_(k) being the angular frequencies of interest; and     -   for broad band noise sources, provision is made to minimize the         following functional f₂ relating to the segments of interest:         f ₂ =∫∥B(ω)∥² dω  (8)

In practice, instead of calculating ∥B∥² it is possible to use a sub-optimal method which consists in replacing the denominator S² by 1. It can then be shown, providing a search is being made for solution directions n_(j) that are sufficiently distinct, that solutions can be obtained that are close to those given by the exact method. The denominator, which cancels when at least two sources coincide, serves to eliminate interfering solutions where a plurality of directions are identical.

When the sub-optimal method comprises broad band processing, it comprises minimizing the following functional f₃: f ₃ =∫det(<T _(k) ,T _(l) *>k, l=0 to M)dω  (9) This magnitude, which is expressed as a linear combination of the cross-correlation functions γ_(ij) of the signals s_(i) and s_(j) taken at points which are themselves a linear combination of the delays: $\frac{{< n_{j}},{c_{i} >}}{c}$

These cross-correlation functions are calculated only for delays having the same order of magnitude as the dimensions of the antenna divided by the speed of sound. Calculation can then advantageously be performed in the time domain as compared with calculation that is usually performed in the frequency domain on the basis of Fourier transforms.

The above-described method thus consists:

-   -   starting from an antenna comprising a plurality of sensors (two         to ten, and preferably two to five) in picking up acoustic or         vibratory information;     -   in amplifying and in filtering the received signals so as to         limit their spectrum and match the sensors in phase and in gain;     -   in digitizing the signals; and     -   in minimizing the projection of the vector s onto the direction         orthogonal to the vectors T_(j), using the following two         calculation techniques a) and b):     -   a) + in obtaining the Fourier transforms of the signals s_(i);         -   + in calculating one of the functionals f₁, f₂, using the             above-defined expressions for the Gramm matrix determinants;             and         -   + minimizing one of the functionals f₁, f₂ depending on the             number of sources to be located;     -   b) + in using a simplified algorithm which consists in         minimizing the functional f₃;         -   + in calculating the cross-correlation functions γ_(ij); and         -   + in minimizing a linear combination of the             cross-correlation functions, depending on the number of             sources.

Once the minimization operation has been performed, the directions n_(j) of the noise sources are determined. Advantageously, it is also possible to recover the characteristics of the noise sources X_(j).

If N=M, i.e. if there are as many sensors as sources, then the system (5) can in general be inverted.

If N≧M, the problem can be reduced to a square system by premultiplying by: ^(t)T* i.e. by the conjugate transposed matrix of T. System (5) then becomes: ^(t) T*.s(ω)=^(t) T*.T.S(ω) I.e. S(ω)=(^(t) T*.T)⁻¹.^(t) T*.s(ω)  (10)

From equation (10), the signals S_(j) can be calculated so as to discover the characteristics of the sources X_(j).

The description below gives an implementation for detecting one noise source (M=1) using N sensors.

This provides: ${{B(\omega)}}^{2} = \frac{\begin{matrix} {{s(\omega)}}^{2} & {\sum\limits_{i = 1}^{N}\;{{\hat{s}}_{i}*(\omega){\mathbb{e}}^{{+ J}\;\omega\frac{{< n},{c_{i} >}}{c}}}} \\ {\sum\limits_{i = 1}^{N}\;{{\hat{s}}_{i}*(\omega){\mathbb{e}}^{{- J}\;\omega\frac{{< n},{c_{i} >}}{c}}}} & N \end{matrix}}{N}$ i.e.: ${{B(\omega)}}^{2} = {\frac{1}{N}\left\lbrack {{\left( {N - 1} \right) \cdot {{s(\omega)}}^{2}} - {\sum\limits_{\underset{k \neq 1}{k,l}}^{\;}\;{{{\hat{s}}_{k}(\omega)}{{\hat{s}}_{1}(\omega)}{\mathbb{e}}^{{j\omega}\frac{{< n},{{c_{1} - c_{k}} >}}{c}}}}} \right\rbrack}$

For broad band sources, it is thus a question of minimizing: ${\int{{{B(\omega)}}^{2}{\mathbb{d}\omega}}} = {\frac{1}{N}\left\lbrack {{\left( {N - 1} \right){\int{{{s(\omega)}}^{2}{\mathbb{d}\omega}}}} - {\sum\limits_{\underset{k \neq 1}{k,l}}{{Re}\left\{ {\int{{{\hat{s}}_{k}(\omega)}{{\hat{s}}_{1}(\omega)}{\mathbb{e}}^{{j\omega}\frac{{< n},{{c_{1} - c_{k}} >}}{c}}{\mathbb{d}\omega}}} \right\}}}} \right\rbrack}$

I.e. writing the cross-correlations between the measured signals as γ_(kl): ${\int{{{B(\omega)}}^{2}{\mathbb{d}\omega}}} = {\frac{1}{N}\left\lbrack {{\left( {N - 1} \right){\int{{{s(\omega)}}^{2}{\mathbb{d}t}}}} - {\sum\limits_{\underset{k \neq 1}{k,l}}\gamma_{k1}^{(\frac{{< n},{{c_{1} - c_{k}} >}}{c})}}} \right\rbrack}$

The problem thus reduces to maximizing: $\sum\limits_{\underset{k \neq 1}{k,l}}{\gamma_{k1}\left( \frac{{< n},{{c_{1} - c_{k}} >}}{c} \right)}$ by varying n which depends on only Np parameters, where Np is equal to the dimension of the space minus 1 unit.

Without further elaboration, it is believed that one skilled in the art can, using the preceding description, utilize the present invention to its fullest extent. The preceding preferred specific embodiments are, therefore, to be construed as merely illustrative, and not limitative of the remainder of the disclosure in any way whatsoever. Also, any preceding examples can be repeated with similar success by substituting the generically or specifically described reactants and/or operating conditions of this invention for those used in such examples.

Throughout the specification and claims, all temperatures are set forth uncorrected in degrees Celsius and, all parts and percentages are by weight, unless otherwise indicated.

The entire disclosures of all applications, patents and publications, cited herein are incorporated by reference herein.

From the foregoing description, one skilled in the art can easily ascertain the essential characteristics of this invention and, without departing from the spirit and scope thereof, can make various changes and modifications of the invention to adapt it to various usages and conditions. 

1. A method of detecting and locating noise sources each emitting a respective signal S_(j) with j=1 to M, detection being provided by means of acoustic wave or vibration sensors each delivering a respective time-varying electrical signal s_(i) with i lying in the range 1 to N, the method comprising: taking the time-varying electrical signals delivered by the sensors, each signal s_(i)(t) delivered by a sensor being the sum of the signals S_(j) emitted by the noise sources; amplifying and filtering the time-varying electrical signals as taken; digitizing the electrical signals; calculating the functional f, such that: ${f\left( {n_{l},\ldots\;,n_{j},\ldots\;,n_{M}} \right)} = \frac{\det\left( {{< {T_{k}(\omega)}},{{T_{l}^{*}(\omega)} > k},{l = {0\mspace{14mu}{to}\mspace{14mu} M}}} \right)}{\det\left( {{< {T_{k}(\omega)}},{{T_{l}^{*}(\omega)} > k},{l = {1\mspace{14mu}{to}\mspace{14mu} M}}} \right)}$  with $\left( {T_{k}(\omega)} \right)_{i} = {\mathbb{e}}^{J\;\omega\frac{{< n_{k}},{c_{i} >}}{c}}$ <.,.> being the scalar product; c_(i) being the vector constructed between the center of gravity of the sensors and the position of sensor i; n_(j) being the unit vector in the direction defined by the center of gravity of the senors and source j; with T₀=s; and with c=the speed of sound; and in minimizing the functional f relative to the vectors n_(j) for j=1 to M in such a manner as to determine the directions n_(j) of the noise sources, wherein ω is angular frequency, <T_(k)(ω), T_(i)*(ω)> is the scalar product between T_(k)(ω), and T_(l)*(ω), J corresponds to the imaginary number in mathematics.
 2. A method according to claim 1, wherein, in order to minimize the functional f when the noise sources are narrow band sources, the method comprises: calculating the Fourier transforms of the signals s_(i)(t) delivered by the sensors; using the expressions for the determinants of the matrices of general term: <T _(k)(ω),T _(l)*(ω)>  to calculate the functional: $f_{1} = {\sum\limits_{k}{{B(\omega)}_{k}}^{2}}$ and after selecting a determined number of noise sources, minimizing the functional f₁ to determine the directions n_(j) of the selected noise sources, wherein B is the noise vector which depends on ω.
 3. A detection method according to claim 1, wherein, in order to minimize the functional f when the noise sources are broad band, the method comprises: calculating the Fourier transforms of the signals s_(i)(t) delivered by the sensors; using the expressions of the determinants of the matrices of general term: <T _(k)(ω),T _(l)*(ω)>  to calculate the functional: f ₂ =∫∥B(ω)∥² dω and after selecting a determined number of noise sources, minimizing the functional f₂ to determine the directions n_(j) of the selected noise sources, wherein dω is the derived term in the integral mathematics formulation.
 4. A detection method according to claim 1, wherein, in order to minimize the functional f, the method comprises: simplifying the expression for the functional ƒ to minimize the following functional f₃: f ₃ =∫det(<T _(k) ,T _(l) *>k, l=0 to M)dω calculating cross-correlation functions γ_(ij) of the signals s_(i) and s_(j); and after selecting a determined number of noise sources, minimizing the functional f₃, wherein <T_(k), T_(l)*> is the scalar product between T_(k) and T_(i)*.
 5. A detection method according to claim 1, wherein, after the minimization operation, the method comprises calculating the source vector: S(ω)=(^(t) T*.T)⁻¹.^(t) T*.s(ω) in order to discover the characteristics of the noise sources, wherein ^(t)T* is the conjugate transposed matrix of T and T is a matrix, S(ω) is the vector of S which depends on ω. 